Photonic crystal optical waveguides for on-chip Bose-Einstein condensates 
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We propose an on-chip optical waveguide for Bose-Einstein condensates based on the evanescent 
hght fields created by surface states of a photonic crystal. It is shown that the modal properties of 
these surface states can be tailored to confine the condensate at distances from the chip surface sig- 
nificantly longer that those that can be reached by using conventional index-contrast guidance. We 
numerically demonstrate that by index-guiding the surface states through two parallel waveguides, 
the atomic cloud can be confined in a two-dimensional trap at about l/im above the structure using 
a power of O.lmW. 

PACS numbers: 42.70.Qs, 42.82.Et, 03.75.-b 



Magnetic microtraps able to guide Bose-Einstein con- 
densates (BECs) have been proposed as a promising way 
to implement on-chip matter- wave waveguides . How- 
ever, this technology is restricted by a fundamental lim- 
itation: the trapping mechanism works for only a sub- 
set of hyperfine states for any given trapping magnetic 
field. It has been shown that this limitation can be cir- 
cumvented in a purely optical trap, where atoms in all 
hyperfine states can be confined . When working with 
an optical trap, one has the possibility of modifying the 
atomic interaction of a BEC by simply varying the ex- 
ternal magnetic field This could be of paramount 
importance for the development of efficient Sagnac in- 
terferometers based on two counterpropagating matter 
waves. Due to the mentioned interest, different kinds of 
dielectric structures have been proposed as optical traps 
SBIilllilSIISIIlIIlIIl However, the extension 
of these ideas to create optical on-chip BEC waveguides 
and circuits is limited. 

One of the simplest schemes one could imagine for on- 
chip waveguiding of BECs would be to use the evanes- 
cent electric field produced by illuminating a high-index 
waveguide [ewav] deposited on top of a uniform dielectric 

substrate of dielectric constant Csufis (being eu,at) > ^subs)- 
In order to have index-guiding in the waveguide, without 
leakage into the substrate, the wave vector k of the guided 
mode has to satisiy: k > (uj / c) ^/e^ubs . This in turn 
means that the decay constant of the evanescent tails into 
air is given by a = ^/W^^uPjc? > {21: / XAiR)\/esubs - 1, 
where Xair is the wavelength of light in air. Because 
Csubs cannot be close to unity (it is typically larger than 
2.0), this is a fast decay that would prevent the balance 
between the repulsion from the surface and the Casimir- 
Polder attraction. For similar reasons, same constraints 
apply for any index-guided scheme that involves a solid- 
state substrate for which Csubs — 1 ~ 1. On the other 
hand, designs without a substrate would have thickness 
< \air, rendering them very fragile and difficult to work 
with. 

In this letter, we move a step forward and we show that 
the use of a photonic crystal (PhC) as a substrate would 
permit to have decay lengths long enough to circumvent 



FIG. 1: Schematic picture of the structure proposed for BEC 
waveguiding. The system consists of a one-dimensional pho- 
tonic crystal with a defect layer on its top. In the defect 
layer two parallel index-guiding waveguides are included. As 
it is sketched in the figure, for a suitable design of the sys- 
tem, the BEC cloud can be confined in the region between 
these waveguides. The reference system used throughout the 
text, together with the geometrical parameters defining the 
structure are also shown. 



the limits imposed by the Casimir-Polder attraction. We 
propose a novel design for optical guiding of BECs based 
on the evanescent light fields created by surface states of 
one-dimensional photonic crystals 15J (ID-PhCs). Sur- 
face states are prohibited from propagating into the PhC 
bulk by photonic bandgap effects, and are prevented from 
leaking into air because of total internal reflection (TIR) 
effects; thus, they are confined to propagate along the 
surface. Surface states of this kind are unique to PhCs. 

Figure 1 shows a schematic picture of the structure an- 
alyzed in this work for on-chip BEC waveguiding. The 
main idea proposed here to get a confining potential both 
in the x- and y-dircctions (see the reference system in Fig. 
1) is the following. First, we create a ID-PhC by a pe- 
riodic sequence of a high dielectric constant {chi) and a 
low dielectric constant (e/o) layers (plotted as red and 
yellow volumes in Fig. 1, respectively). The thickness 
of each layer is given by Whi and Wio- In addition, we 
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FIG. 2: (a) Projected band structure of the semi-infinite one- 
dimensional photonic crystal (PhC) considered in this work 
[see the profile of dielectric constant in panel (b)] . Red shaded 
areas show the regions where there exist extended states inside 
the PhC. Grey shaded area corresponds to eigenmodes that 
are propagating in air. Orange, blue and magenta lines rep- 
resent the dispersion relations of the surface states for several 
values of the thickness of the top layer. Green shaded region 
corresponds to possible confined states of any defect layer if 
it lies on top of a uniform dielectric substrate with e = ejo. 
(b) Dielectric constant (black line) and electric field inten- 
sity profiles along the direction perpendicular to the interface 
PhC-air. Three surfaces states with frequencies close to the 
light line are displayed for Wd— 0.07a. 



define a as the periodicity of the PhC. Secondly, we intro- 
duce a defect in this structure by reducing the thickness 
of its top layer to Wd. The electric field of this defect 
mode is confined both in y> and y< regions, i.e., 
this mode is a surface state. By coupling into the de- 
fect layer a laser beam whose frequency is blue-detuned 
from a certain transition of the atoms forming the con- 
sidered BEC, the condensate will experience a repulsive 
potential from the surface. This potential combined with 
the Casimir-Polder interaction and the force exerted by 
the gravity (which pushes the atoms towards the surface) 
will trap the BEC in the minimum of the potential in the 
y-direction. 

Once a suitable confinement in the direction normal 
to the air-PhC interface is obtained, the next step is to 



design a practical mechanism to confine the BEC in the 
X direction. We have found that one convenient way to 
achieve this is to create two index-guided parallel waveg- 
uides for these surfaces states, (represented as blue vol- 
umes in Fig. 1, where Ww and Hy^, stand for the width 
of the waveguides and the separation between them, re- 
spectively). 

In order to optimize the potential created by the struc- 
ture described in the previous paragraphs we first analyze 
the case in which the PhC is uniform in the x-direction 
(i.e. without waveguides). This is a good starting point 
because the decay of the electric field in the y-direction 
will be almost the same as in the system with the two 
waveguides 1 1 6j. 

Red regions in Fig. 2(a) show the projected bands of 
a semi-infinite ID-PhC computed by means of t he p lane- 
wave expansion method to Maxwell's equations [l7|. We 
have taken the following values for the parameters defin- 
ing the structure e?ii=12.25, e;o=2.25 and Wio=0.7a, as 
they are typical values that could be used in the exper- 
iments. In addition we have assumed e=l in the region 
above the PhC. The grey shaded region represents modes 
that can propagate in air. In order to show how the dis- 
persion relation for the surface state changes with the 
thickness of the defect layer, we have also plotted in Fig. 
2(a) the results for three different values of Wd (ma- 
genta, blue and orange lines correspond to 1^^=0.040, 
Wd=0.07a and Wd=0.1a). Any of these values for Wd 
could be chosen for our design; for definiteness, from now 
on we take Wd=0.07a. 

Now, let's turn our attention to study how the elec- 
tric field profile depends on the frequency of the de- 
fect mode. Figure 2(b) shows a cross-section in the y- 
direction of electric field intensities of the surface states 
corresponding to three representative frequencies near 
the light line. Blue, green and red lines render the results 
for wa/27rc=0.214, wa/27rc=0.225 and uja/2nc=0.235, re- 
spectively. The profile of the dielectric constant along the 
y-direction of the structure has also been plotted. As it 
is expected, the decay length in air increases as we get 
closer to the light hne (i.e. as the frequency decreases). 

At this point, it is worth taking a closer look at the 
advantages of the type of systems studied in this work 
compared with other kinds of rigid dielectric structures 
that are not based on PhCs. As it is shown in Fig. 2(a), 
the dispersion relation of the defect mode of a PhC allows 
us to choose a frequency as close to the light line as we 
need in our atomic trap design; thereby, the decay length 
in air can be tuned to be arbitrarily long. However, this 
property could not be obtained by means of any TIR- 
based waveguide placed on top of a dielectric substrate, 
as in that case the frequencies of the confined states ujd 
are such that uJd < ck/^/e^, where is the dielectric 
constant of the substrate. This point is illustrated in 
Fig. 2(a), where the green shaded area corresponds to 
the region where confined states can exist for the case 
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Eg = eio- As we can see in that figure, all of those states 
are far from the light line, so their decay length in air is 
short. 

Now, before proceeding with the analysis of the full 
waveguiding structure, let's examine the confining atomic 
potential created by the ID-PhC with a defect described 
previously, but without waveguides. If we illuminate this 
defect layer with a laser beam whose frequency u is blue- 
detuned with respect to a certain atomic transition of fre- 
quency Wo, then an atom located above the PhC will ex- 
perience a potential given by the sum of the optical-dipole 
{Uop), gravitational {Ugv), and Casimir- Polder interac- 
tions (Ucp), that is, J7(r) = J7op(r) + Ugy{r) + UlpI^ 
The optical dipole contribution can be written as jl4lll8| 



hT^ I{r) 



8A L, 



(1) 



where Ps and Igat are two parameters whose value de- 
pend on the atomic transition we are considering, F is 
the corresponding linewidth and A = uj — ujo (note that 
in deriving Eq.(l), it has been assumed that |A| << cjq)- 
/(r) = ieo|E(r)p is the electric field intensity. Since 
we are currently assuming that our structure has trans- 
lational symmetry both in the x and z directions, we 
can write J{r) = Iq exp(— 2y/A), where A defines the 
evanescent decay length of the electric field and Iq is the 
intensity just above the illuminated layer. 

The second term in Eq.(l) represents the gravitational 
potential and its given by Ugv{r) = mgy, where m and 
g stand for the mass of each atom and the gravita- 
tional acceleration, respectively. Finally, the last term 
of Eq.(l) corresponds to the so-called Casimir-Polder po- 
tential Ucp{v) = —C4/y^, where the constant C4 is given 
by in 
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In the above expression, a is the atomic polarizability 
and the value ip{etii)=0-G9 has been computed using the 
analytical expressions for the function ip{ehi) given by 
Yan et al. fig, [H] 

In order to show the feasibility of the experimental 
realization of the system presented in this work, we con- 
sider the case of Cs atoms, since the confinement of a 
condensate of Cs atoms in a similar structure has been 
demonstrated 0. As we are considering the D2 res- 
onant transition, we must take the following values in 
Eq.(l), ps=2/3, /^at=l.lniW/cm2, wo=27rx3.5x lO^Hz 
(Ao=857.1nm), and F=27rx5.3x 10^ Hz. Regarding the 
external illumination, we assume that it is produced by a 
Tirsapphire laser at a;=27rx3.6x lO^'^Hz. In addition, to 
get an optical potential with a similar depth to that ob- 
tained in Ref. 13, we suppose that the maximum electric 
field intensity in the air region very close to the PhC-air 
interface is /o=2. 67x10^ mW/cm^. 



TO 

O 

O 
Q. 

O 

E 
o 

< 




A=0.6nm 
A=0.3nm 
A=0.2nm 
A=0.1nm 



1 2 3 4 5 6 
Distance from the surface (|am) 



FIG. 3: Main panel: Total atomic potential for the case of Cs 
atoms as a function of the distance to the dielectric surface. 
The results for several decay lengths of the evanescent light 
field are shown. The optical dipole potential is assumed to be 
created by a Ti:sapphire laser at a wavelength of 839nm. All 
the cases plotted in this figure were computed for the same in- 
tensity 7o=2.67xlO''mW/cm^ evaluated just above the PhC- 
air interface. Inset: Total, optical, gravitational and Casimir- 
Polder potentials (blue, green, red and black, respectively) for 
the case of A=0.3/im. The minimum for the total potential is 
found at i/o=l-l/j,m. 



Figure 3 plots the total atomic potential as a function 
of the distance from the surface for several values of the 
decay length A ranging from 0.1 to 0.6/xm. As it is shown 
in this figure, a decrease in A leads to a shift of the min- 
imum of the potential towards the surface. However as 
can be seen, there exists a minimum distance from the 
surface where the condensate can be trapped, since if we 
take A <0.2/xm the total potential does not have a min- 
imum. This limit is established by the fact that in the 
region where the repulsive optical potential is not negligi- 
ble for A <0.2/im, the Casimir-Polder attractive interac- 
tion dominates, preventing the formation of a minimum 
in the total potential. In the case of the corresponding 
TIR waveguide design with a substrate with e =2.25 one 
would obtain a maximum decay length A <0.12yLtm; thus, 
given the rapid variation of Casimir-Polder interaction as 
a function of distance, this demonstrates explicitly that 
a TIR scheme can not be made to work for the current 
purpose, and that PhCs indeed have to be used. 

Since we are interested in bringing the BEC as close 
as possible to the surface, we will choose A=0.3/im. This 
value corresponds to a state of wavevector A;a/27r=0.27 
of the defect band, and thereby a value of the periodicity 
a=0.l9fim (see Fig. 2(a)). From now on, we will assume 
this value of k, unless otherwise stated. In that case, the 
minimum of the potential appears at ?/o=l-lA*ni- 

It is also interesting to analyze the contribution of each 
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FIG. 4: 2D total atomic potential produced by the structure 
sketched in Fig. 1 in the xy plane for the case of Cs atoms. 
Color scale codes the potential in units of fiK. The minimum 
is found at (a^o, 2/o) = (3.3,l.l)/im. Left and right insets show 
cross-sections of the potential along the lines yo=l.lij.m and 
xo~S.3fim, respectively. 

term appearing in Eq.(l) to the total atomic potential. 
Blue, green, red, and black lines in the insert of Fig. 3 
show the total, optical, gravitational and Casimir-Polder 
potentials, respectively, for A— 0.3/im. This figure reveals 
that at small distances the Casimir-Poldcr potential gov- 
erns the behavior of the total potential while for distances 
slightly larger that i/q, the total potential is dominated 
by the gravitational interaction. 

Finally, let us focus on how to obtain a suitable confin- 
ing potential in the x-direction, needed to create an op- 
tical BEC waveguide. As we described above, we could 
solve this problem by including in the defect layer two 
index-guided waveguides of dielectric constant e^, with 
> Shi (see schematics in Fig.l). For definiteness, we 
take eu, = 18 The mode we explore is the even mode 
of this system of two coupled waveguides. In order to 
reproduce faithfully the features that could be found in 
the experiments, we have simulated the structure u sing 
the Finite-Difference Time-Domain (FDTD) method |23j. 
The results obtained with this method come from exact 
solutions of the 3D Maxwell equations, except for the 
numerical discretization. 

As we described in the introduction, one of the main 
advantages of our proposal is that there are several 
parameters that can be changed to our convenience. 
Through the previous calculations, we have fixed the val- 
ues of the thickness of the defect layer Wd, as well as the 
wavevector k and the frequency to considered. In addi- 
tion, now for simplicity, we fix the width of the waveg- 
uides to Ww—'^.a. Therefore, the remaining parameters 
in our optimization process are the distance between the 
waveguides and the power Pq carried by the even 



mode they create. We now look for the values of [H^^Pq) 
such as the curvature of the potential at its minimum is 
similar both in the x- and the y-directions. This con- 
dition would be useful if we want to guide a cylindrical- 
shape BEC along our waveguide. In addition, notice that 
we require a low value of Pq in order to have an efficient 
guiding system. 

Following the program described above, we have found 
a suitable configuration for _ffu,=3.2/im and Po=0.1mW. 
Figure 4 shows a contour plot of the corresponding two 
dimensional total atomic potential in the xy plane. The 
color scale codes the potential in units of ijlK. This fig- 
ure shows a minimum at (a;o,j/o) — (3.3, l.l)/im, be- 
ing the center of the waveguides located at xi=1.5/im 
and X2=5.1/J.m. Left and right insets plot cross-sections 
along yQ—l.l^Tci and a;o=3.3/im, respectively. From these 
cross-sections, we find that the trap frequency at the 
minimum along the x-axis is Wa;=1233Hz, while along 
the y-axis it is ujy—QA2Hz. This corresponds to a ratio 
L0x/ojy=l.'d2, which confirms that it is possible to get a 
confining potential with similar curvatures in both trans- 
verse directions using this configuration. Notice that, 
with the same power, we could get other shapes for the 
BEC by changing the ratio uJx/'jJy through the variation 
of H^. 

In conclusion, we have proposed a new class of optical 
BEC waveguides based on illuminating two waveguides 
built into a defect of a ID-PhC. We have demonstrated 
that a power as low as O.lmW is enough to guide a BEC 
of Cs atoms through the suggested waveguide at about 
l//m above the structure. The kind of system proposed 
in this work could be considered as a basic ingredient for 
the development of novel BEC on-chip devices. 
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